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We theoretically demonstrate that, in a ferromagnet/paramagnet bilayer, a magnon instability 
accompanied by a gigahertz microwave emission can be driven simply by means of a temperature 
bias. Employing many-body theory for investigating the effects of a phonon heat current on the 
magnon lifetime, we show that the magnon instability occurs upon the suppression of the Umklapp 
scattering at low temperatures, leading to microwave emission. The present finding provides crucial 
information about the interplay of spin current and heat current. 

PACS numbers: 85.75.-d, 72.25.Mk, 75.30.Ds 


I. INTRODUCTION 

The concept of the spin-transfer torque has enabled 
the establishment of a quantum-mechanical way of con¬ 
trolling spin-based electronic devicesi^— . Although these 
spintronic devices are driven by spin-polarized electric 
currents, a recent trend among other important issues in 
spintronics is the thermal manipulation of these devices. 
Indeed, the spin Seebeck effecl^P— has enabled a thermal 
spin injection from insulating magnets into paramagnets, 
and the discovery has opened up a new research field 
known as spin caloritronics^. The field has continuously 
expanded and given rise to a number of intriguing phe¬ 
nomena emerging from the interplay between spin and 
healjiSrJ^. 

The basic building block of the spin Seebeck device is 
a bilayer composed of a ferromagnet and a paramagnetic 
heavy metal. This simple ferromagnet/paramagnet bi¬ 
layer differs from the conventional spin-transfer torque 
device in that unlike the latter system, the former does 
not possess the so-called spin valve structure and thus the 
device structure can be much simpler. Given the increas¬ 
ing attention focused on the thermal effects in spintronics 
and the versatility of the simple ferromagnet/paramagnet 
bilayer, theoretical clarihcation of the thermal manipula¬ 
tion of ferromagnet /paramagnet bilayersi^^— is of crucial 
importance. 

In this work, we theoretically address the issue of 
thermally driven microwave emission from the ferromag¬ 
net/paramagnet bilayer system. Using a diagrammatic 
perturbation calculation with respect to the external 
temperature bias as well as the interfacial magnetic inter¬ 
action, we show that phonon heat current has remarkable 
influences on the magnon lifetime in the bilayer via the 
phonon drag mechanism, enabling thermal manipulation 
of the Gilbert dampin g. T he theory explains the exper¬ 
iment reported in Ref. [Ij, where thermal control of the 
Gilbert damping in a prototypical Y 3 Fe 50 i 2 /Pt bilayer 
was demonstrated. Moreover, we predict that the Gilbert 
damping constant in the same bilayer system can become 
negative provided that Umklapp scattering of phonons 


is sufficiently suppressed at low temperatures, thereby 
causing an instability in the magnon system, leading to 
the microwave emission. 

The present paper is organized as follows. In Sec. |nl 
we define our theoretical model to deal with a ferromag¬ 
net/paramagnet bilayer system and outline how the en¬ 
hanced Gilbert damping in the system is calculated. In 
Sec. mi we introduce a temperature bias to our model 
and discuss its effects on the enhanced Gilbert damp¬ 
ing. Based on this result, in Sec. El we predict that 
thermally driven microwave emission from the bilayer is 
possible at low temperatures. Finally, in Sec. El we dis¬ 
cuss the physics behind the predicted microwave emission 
and then summarize our result. 


II. MODEL AND ENHANCED GILBERT 
DAMPING 


The system under consideration is a bilayer composed 
of an insulating ferromagnet {F) and a metallic param¬ 
agnet {P) under a temperature bias (Fig. IT]). An external 
static magnetic held iJo is applied in the lateral direction, 
and the anisotropy held iJan is assumed to be in the same 
direction, with its strength being much weaker than Hq. 
The physics of F and P are, respectively, described by 
the localized spins and the spin accumulation, and the 
interaction between them is given by the s-d exchange 
coupling at the interface. Since, in our approach, the 
nonequilibrium dynamics of these quantities play a key 
role, it is convenient to employ a diagrammatic pertur¬ 
bation approac h^^d^ formulated in terms of the quantum 
action. We use the matrix form of the propagators repre¬ 
sented in terms of the retarded, advanced, and Keldysh 

components2£. 
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FIG. 1: (color online) Schematic of the magnetic bilayer 
and a physical picture of thermally driven microwave emis¬ 
sion, where Hq is an applied static magnetic held and Jfan 
is the anisotropy held (/fan -C Hq). (a) In the absence of 
a temperature bias, a precessing magnetization in the ferro- 
magnet F, or a magnon, creates spin current /p^^p pumped 
into the adjacent paramagnet, P. This effect always results 
in an enhancement of the Gilbert damping constant because 
the magnons lose their spin angular momenta, (b) When the 
paramagnet is heated, the pumped spin current tends to dif¬ 
fuse back into the ferromagnet because of the drag effect of 
the phonon heat current Qph, and hence produces a forced 
spin backhow. This spin backhow reduces the Gilbert damp¬ 
ing constant, as the magnons gain additional spin angular 
momenta, thereby causing a magnon instability, leading to 
microwave emission. 


Our starting point is the following action! 

Ao = f ^a]j(t)Gq^(i5t)aq(t) 

JC q 

k 

q.k 

where the integration is performed along the Keldysh 
contour C. The first term in Eq. m describes the dynam¬ 
ics of magnons, represented by the Holstein-Primakoff 
operators Oq and a^. Here the retarded part of the 
magnon propagator is given by G^{uj) = l/{uj — ujq + 
iof^^^o;), with Wq and being the magnon excitation 
energy and the intrinsic Gilbert damping constant, re¬ 
spectively. In the following, the uniform magnon fre¬ 
quency is denoted as ojq. The second term in Eq. (HD 
describes the dynamics of the paramagnetic spin den- 
sity 2 ^ (spin accumulation), represented by = Skiis^, 
and the retarded part of the propagator is given by 
Xk(^) — Xp /(1 + — iwTgf), where xp is the uni¬ 

form paramagnetic susceptibility, A is the spin diffusion 
length, and Tgf is the spin-flip relaxation time^^. The 
third term in Eq. (HD describes the s-d exchange interac¬ 
tion at the F/P interface, where Jsd(k — q) is the Fourier 


transform of Jsd Sroeinterface - To) with Jsd and Vg 
being the strength of the s-d interaction and the cell vol¬ 
ume, respectively. 

We begin by explaining how the damping of magnons 
in this bilayer system, conventionally represented by the 
effective Gilbert damping constant, is calculated in our 
approach. Guided by the experimental observation that 
the use of a poor spin sink such as Cu makes the tem¬ 
perature bias effect invisibl o^^A’^ , we assume that P is a 
good spin sink with sizable spin-orbit scattering, such as 
Pt. Since the attachment of a good spin sink introduces 
a new source of spin dissipation, the Gilbert damping 
constant in the bilayer is given as a sum of the intrinsic 
contribution (defined merely by F), and an enhanced 
Gilbert damping constant i5a caused by the attached P, 
i.e., a = +5a. As discussed in Ref. 1^. the enhanced 

Gilbert damping constant is in general calculated as 
5a = —(l/wo)ImE^(wo), where T,q{ujq) is the selfenergy 
for the uniform-mode magnon that is defined by the fol¬ 
lowing Dyson equation^!: l/g^{uj) = 1 /Gq{uj) — S^(w) 
with gg{oj) being the renormalized propagator of the 
uniform-mode magnon. 

The enhanced Gilbert damping constant 5a is calcu¬ 
lated from the magnon selfenergy shown in Fig. HJa). 
Performing a perturbative approach with respect to Jsd, 
we obtain 


A(2) 

k ° 

where ((J^d)) = ‘2‘J‘idSgNipt/{NpNF), with h, Wnt, Np, 
and Np being, respectively, the Planck constant divided 
by 27r, the number of localized spins at the interface, and 
the number of lattice sites in F and P. Such enhance¬ 
ment of the Gilbert damping has been confirmed experi¬ 
mentally by linewidth measurements using ferromagnetic 
resonance (FMR) ^^i^^ . 

The result can be interpreted from the viewpoint 
of spin transfer across the interface^li^, which is now 
termed “spin pumping”. Because the exchange inter¬ 
action at the interface conserves the total spin, the en¬ 
hanced Gilbert damping, or spin angular momentum loss, 
should be accompanied by a spin transfer from F into P, 
where additional spin d issip ation occurs. Following the 
procedure given in Ref. l2lL a spin current, that 

is pumped from F into P under the FMR condition is 
shown to be intimately connected to i5a through the re¬ 
lation 


/P"“P = ,5a 


(a(o))2wo 


( 3 ) 


where 7 and are, respectively, the gyromagnetic ra¬ 
tio and the amplitude of the oscillating microwave field. 
Since a magnon carries a spin angular momentum —h, 
this effect pumps negative spin angular momenta into P 
[Fig-UDa)]. The existence of the pumped spin current has 
been demonstrated experimentally using the inverse spin 
Hall effect^. 
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The perturbation approach with respect to the s-d 
interaction is justified when Jsd is of the order of sev¬ 
eral tens meV such as Y3Fe50i2/Pt systems^i^. In¬ 
deed, we find that the expansion parameter character¬ 
izing the strength of the higher-order terms compared 
to the lowest-order term is given by e = J^^xp/J ex, 
where Jex is the exchange integral of F. Because our 
definition of spin susceptibility differs from the usual 
magnetic susceptibility by a factor (/is: the Bohr 
magneton), the expansion parameter is calculated to be 
Xp ^ ^stonerZ-^F, where Astoner and Ep are respec¬ 
tively the dimensionless Stoner enhancement factor and 
Fermi energy of P. Therefore we have an estimate e ^ 
(Jsd/Tex)(^Stoner>/sd/-EF) < 1 , where we have used the 
fact that, in the case under consideration, Jsd/Jex ~ 1 
whereas ^stonerTsdZ-E'F -C 1 even if we consider a mod¬ 
erate strength of the Stoner enhancement factor of Pt, 
^Stoner ^ 5 ^. This estimate justifies the perturbation 
approach with respect to Jgd in Y3Fe50i2/Pt systems. 


III. EFFECTS OF TEMPERATURE BIAS 

Now, we discuss the effects of a temperature bias on 
the enhanced Gilbert damping. We consider a situation 
where F and P are in contact with heat baths with tem¬ 
peratures Tcoid and Thot, and where the thermal resistiv¬ 
ity and the thickness of P are considerably smaller than 
those of F, such that the temperature profile of the bi¬ 
layer behaves as shown in Fig.[ 2 ])b). It is well known that 
the temperature is a statistical property of the system; 
therefore, no Hamiltonian exists to describe the effects 
of a temperature bias^^iM. To consider the thermal ef¬ 
fects in our approach, a model is introduced in which the 
system is divided into three temperature domains with 
local equilibrium temperatures of Ti = (Thot + 3 Tcoid)/ 4 , 
T2 = ( 3 rhot + Tcoid)/4, and T3 = Thot- We assume that 
each domain is initially in local thermal equilibrium with¬ 
out interaction with the neighboring domains, and then 
switch on the interactions between them, and calculate 
the nonequilibrium dynamics of the system. This method 
has been proven to be quite efficient for a theoretical de¬ 
scription of the spin Seebeck effect^. 

The size of the temperature domains may be given 
by a length scale defining local equilibrium temperature, 
which is set by the mean free path of acoustic phonons^. 
In the case of a widely used Y3Fe50i2, if we neglect 
so-called subthermal effects by which low-energy long- 
ranged phonons carry most of the heat, it is of the order 
of 1 nm at 300 K and 1 ^m at 50 K^. If we take into 
account the subthermal effects, the length is expected to 
be much longer. Apart from this estimate, however, the 
important point of the following calculation is that the 
result scales with the temperature gradient, instead of 
the temperature difference itself [see the paragraph be¬ 
low Eq. m]- Thus, the following result is independent 
of the choice of the domain size. 

We apply this technique to the present problem and 



FIG. 2: (color online) (a) Diagrammatic representation of 
the basic process that gives the enhanced Gilbert damping, 
which is free from the temperature bias effect. Here the bold 
solid line and the double solid lines with a pair of opposite ar¬ 
rows, respectively, represent the propagators of the magnon 
and the spin accumulation, (b) Temperature profile of the 
system, where the paramagnet P is in contact with a heat 
source of temperature Thot and the ferromagnet F is in con¬ 
tact with a heat sink of temperature Tcoia. In our model¬ 
ing, the continuous temperature profile (thick solid line) is 
replaced by three discrete local temperatures of Ti, T 2 and 
T 3 . (c) The phonon-drag process that gives the thermal con¬ 
trol of the Gilbert damping. Here the dashed line represents 
the phonon propagator, (d) The Feynman diagram for the 
longitudinal spin Seebeck effect discussed in Ref. [h] . 


calculate the enhanced Gilbert damping constant. We 
then find that, as long as model m consists merely of 
magnons and spin accumulation, there is no tempera¬ 
ture bias effect on the enhanced Gilbert damping, even 
though there is a finite temperature jump T3 — T2 7^ 0 at 
the F/P interface in our model calculation. That is, even 
if we evaluate the diagram in Fig. [ 21 [a) under a temper¬ 
ature bias (Thot Tcoid), no effects arise from this tem¬ 
perature bias. This result originates from the fact that 
the enhanced Gilbert damping is a purely “mechanical” 
process^ that is not represented by magnon distribution 
function. For the thermal bias to play a role, we need 
to consider a process that involves, in addition to the 
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magnons and spin accumulation, a third degree of free¬ 
dom that carries heat and interacts with the magnons or 
spin accumulation. 

To consider the temperature bias effects, we thus focus 
on acoustic phonons as the third degree of freedom. In a 
widely used Y 3 Fe 50 i 2 /Pt bilayer system, it is well estab¬ 
lished that acoustic phonons give a quite sharp peak to 
the spin Seebeck effect, owing to the phonon-dra g rn ech- 
anisnt^^— at low temperature. Following Ref. |^, we 
consider the following action: 


4lph 


k 


+ XI '^K-BK(t)Sk+K W}, 

k,K 


(4) 


where B^ ~ 4" is the creation operator for 

phonons, and D^{v) = (v — i/rph)“^ — {v + vj<i + 
i/'T’h)”^ is the retarded part of the phonon propagator, 
with Tph and r'K = Vp\iK being the phonon lifetime and 
energy; here the phonon propagator depends on the un¬ 
derlying material. The first term describes the dynamics 
of the acoustic phonon, and the second term describes 
the interaction between the phonon and the spin ac¬ 
cumulation in P, where the explicit expression of the 

vertex is given by Tk = (Ss-p/fi) with 

Mion and g^-p being the ion mass and the strength of 
the spin accumulation-phonon interaction^, respectively. 
Note that in Eq. ®, effects of the phonon-phonon Umk- 
lapp scattering is phenomenologically taken into account 
through the temperature dependence of the phonon life¬ 
time Tph, which is know to be largely suppressed at room 
temperature by Umklapp processes. Besides, the spin 
accumulation-phonon Umklapp scattering is ill-defined 
because the spin accumulation is a hydrodynamic dif¬ 
fusive mode. Therefore, Eq. o is valid even at room 
temperature. 


In the following, we show that, in the presence of 
phonon heat current concomitant to the temperature 
bias, the total Gilbert damping a can be written as 


a = -I- (5a -I- Sa'AT, (5) 


where the last term represents the thermally driven 
change in the Gilbert damping. Notably, as shown be¬ 
low, the term Sa'AT is proportional to the phonon heat 
current, Qph, a situation that is accounted for only by 
nonequilibrium phonon dynamics. 

The dominant contribution to 6a' AT at low tempera¬ 
ture comes from the process shown in Fig. [UJc) , where the 
acoustic phonons carry the heat current Qph and drag 
the spin accumulation. As before, we have a relation 
Sa'AT = — (l/wo)ImSQ^(a;o) for the present problem. 


where 




j2 q 

V2N,N} A Tk,[xuM1 


X i / Xk-K3(w-J^)5T>^3M 
J y 


( 6 ) 


is the magnon selfenergy corresponding to the diagram 
in Fig. He), and SD^^ {y) is given by 


, hi/ s \ 


(7) 


We evaluate the integral over 1 / by picking up the domi¬ 
nant condribution proportional to Tph, and obtain 




72 o /\r(2i3) 


iTphi X 

k.Ki.Ka.Ks 

X [Xk_(‘^-)+Xk+(w+)]Ini£>g3(7/K2)] 


Nl,Nl 


tL[x^M]^ 




'2kBT3' 


'2kBTi 


( 8 ) 


where k± = k±Ki, uj± = ujAiz-k^, and we have used the 
relation 1tbD^{—v) = Then, after a tedious 

but straightfoward calculation, we have 


SaAT - —SaxpT^sLrp^ j 


Kh(^F))'(fcBT)MAT7 


C 55 

(9) 


where i/p is the Debye frequency of E, 


f'hi'p / k'^T 


C 5 = 


dx 


4sinh^(a;/2) ’ 


( 10 ) 


and the quantity L is defined below Eq. (IA2I1 in Ap¬ 
pendix 1^ 

Note that the above result [Eq. (H] scales with the 
temperature gradient instead of the temperature differ¬ 
ence itself, which can be seen from the fact that Sa and 
L in Eq. m are, respectively, proportional to 1 /Np and 
■ Because A^f = L^LyL^/a^ and = LyLz/o? 

where a is the lattice spacing and LxLyLz is the volume 
of F with the a;-axis chosen along the temperature gra¬ 
dient, Eq. (H is proportional to a{AT/L^) = aVT. 

The right hand side of Eq. (H can be rewritten in terms 
of phonon heat current Qph across the F/P interface. 
Using Eqs. dMI and (IA5I1 in Appendix I3 the thermal 
contribution to the enhanced Gilbert damping can be 
expressed as 


Sa'AT = ^xl^lsBSa Qph, (11) 

TT^ 
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FIG. 3: (color online) The temperature dependence of the 
total Gilbert damping constant a = + Sa + Sa'AT calcu¬ 

lated for a Y 3 FesOi 2 /Pt bilayer under a temperature bias of 
AT = 3 K. Here the non-thermal terms and Sa are taken 
from Ref. 0. while the thermally driven term 5a'AT is cal¬ 
culated using Eq. Hi). Inset: The temperature dependence 
of the spin Seebeck signal S'sse taken from Ref. [39l . which is 
used to calculate 5a' AT via Eq. (nni). 


where ^ Ek ^'K/(i^K + l/^ph) is 

the effective coupling constant between the phonon and 
the spin accumulation, and a dimensionless coefficient B 
is given by 

^ ^ Cs/Cg 

{hvpYPphiPp){plY,iPp)Y (^ 3^)3 

Equation (El) is the main result of the present paper. 
As introduced before, the quantity Qph = —ATphAT in 
Eq. dill) represents the phonon heat current under the 
temperature bias AT = Thot — Tcdd, where the expres¬ 
sion for ATph (> 0) is given by Eq. dA5l) . Note that, 
within a linear response with respect to AT, the phonon- 
drag process discussed here gives the leading contribution 
to Sa' at low temperature, which also well explains the 
phonon-drag peak in the thermal conductivity and the 
spin Seebeck effect in LaY 2 Ee 50 i 2 ^. 


IV. MICROWAVE EMISSION DRIVEN BY 
TEMPERATURE BIAS 

The most remarkable prediction of the present the¬ 
ory is that the total Gilbert damping constant a in 
Eq. ([5]) can become negative, signifying an instability 
in the magnon system and giving rise to gigahertz mi¬ 
crowave emission. The prediction is made as follows. 
Eirst, according to Eq. ED, the sign of Sa'AT can be 
negative when P is hotter because of the condition that 


Qph = -A'ph(Thot - Ecoid) < 0. Second, the magni¬ 
tude of Sa'AT is significantly enhanced at low tempera¬ 
tures where the phonon thermal conductivity ATph has a 
phonon-drag peak^, thus compensating for the positive 
contribution from -I- Sa and leading to a negative 
Gilbert damping constant. 

To judge the feasibility of the predicted microwave 
emission, we calculate the temperature dependence of a 
using information given in the literature. In doing so, we 
first limit ourselves to a typical Y 3 Ee 50 i 2 /Pt bilayer sys- 
teroi^^— , in which the reduction of the Gilbert damping 
constant under a temperature bias is actually observed. 
We next use the fact that the longitudinal spin Seebeck 
effect in the same system is governed by a quite similar 
phonon-drag process^EilL, shown in Eig. HKd). Use of the 
longitudinal spin Seebeck signal allows us to eliminate ex¬ 
perimentally unknown parameter YpT^jj from Eq. dill) . 

We can calculate the spin current generated by the 
longitudinal spin Seebeck effect, /|®®, following the pro¬ 
cedure given in Ref. 10 (see Appendix iBl for details). The 
result is 


rSSE 

-'S 


= SaN^47r^LTp^r!^^^A^^ 

/AT 


{hvpY 


\ T 


C 3/2 C6, (13) 


where Pmag(w) is the density of states of magnon, wp = 
27r^ Jex'S'o /h with Jex being the exchange energy of T, and 
C 3/2 and cq are respectively defined by 




C 3/2 


dx 


, 3/2 


sinh(a;/2) ’ 


nhiyp jk^T 


C6 = 


dx- 


(14) 

(15) 


Jo 4sinh(a::/2) 

By comparing Eq. ([5]) with Eq. m, we finally obtain 
Sa'{T)AT = (l/a;o)C(T)/|SE(T), (16) 


where a dimensionless coefficient C{T) is given by 

C{T) = 3 C5 

7rpmag(wD)(wDTsf)2 C3/2Co'^ 


Note that the above result is derived within the linear- 
response calculation and thus the result is always linear 
with respect to the temperature bias AT. 

An important consequence of Eq. is that it gives 
the relation 


Sa'{T)AT 




( \ 


(18) 


where Tr is the room temperature. From this, we can es¬ 
timate Sa' {T)AT, because the temperature dependence 
of the ratio i^|®®(T)//|®®(TR) is given in Ref. HI, the 
value of Sa'{T^ijAT is known from Ref. M and the value 
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of C{T)/C{Tji) can be calculated numerically. The va¬ 
lidity of this procedure is firstly checked by conmaring it 
with an experimental result reported in Ref. Il4l where 
a 17% reduction in the Gilbert damping constant in a 
Y 3 Fe 50 i 2 /Pt bilayer is observed at room temperature 
when the Pt side is hotter by an amount AT = 20K (see 
Appendix ICl). 

Figure [3] shows the temperature dependence of the 
total Gilbert damping constant a, calculated for a 
Y 3 Fe 50 i 2 /Pt hybrid system, where the non-thermal 
terms and Sa are taken from Ref. [13 and assumed 
to be temperature-independent. Glearly, we find that a 
negative Gilbert damping constant is achieved at around 
T = 50 K under a temperature bias of AT = 3 K, which 
signifies an instability in the uniform-mode magnon dy¬ 
namics. 

This magnon instability with the negative total Gilbert 
damping constant indicates that the energy gain due to 
the work done by the temperature bias exceeds the en¬ 
ergy loss due to the non-thermal Gilbert damping con¬ 
stant, -I- (5a, such that the excess energy dissipates 
in the form of emitted microwave radiation. This is sim¬ 
ilar to the case of microwave emission from spin valves 
that are driven by a spin-transfer torqu o'^^d^ . Unlike 
spin valve systems, however, magnetization reversal does 
not occur; for the magnetization reversal to be real¬ 
ized, two well-defined energy minima in the free energy— 
corresponding to parallel and anti-parallel states in the 
case of the spin valve—are required. In the present case, 
on the other hand, because we are dealing with a situ¬ 
ation where the external magnetic field is considerably 
larger than the in-plane anisotropy field as stated in the 
beginning of this section, the system possesses only one 
energy minimum, which corresponds to the direction of 
the external magnetic field. In such a high-held region, 
the magnon instability signihes only microwave emission, 
even in the spin valve system^. 

Therefore, this calculation shows that we can drive a 
gigahertz microwave emission with work done by an ap¬ 
plied temperature bias. Because the uniform mode has 
the narrowest linewidth in the spin-wave resonance spec¬ 
trum in the Y 3 Fe 50 i 2 /Pt system^, the instability is ex¬ 
pected to originate from uniform precession. However, 
since a magnetostatic surface mode may have stronger 
coupling to the Pt layer because of its surface-localized 
nature, an instability of this mode prior to the uniform 
precession may also be possible. Although the microwave 
emission in an actual experiment is expected to be con¬ 
comitant with the dynamically precessing state involving 
highly nonlinear processes, the present theory deals only 
with linear stability analysis; the precise description of 
such a dynamically precessing state is beyond our scope. 

It is important to note that the predicted microwave 
emission can be separated from background thermal ra¬ 
diations by measuring the power spectra of emitted mi- 
crowaves^Si^. The predicted microwave emission has 
peaks only at several gigahertz frequencies, which is spe¬ 
cific to magnon dynamics (see Fig. 3 in Ref. By 


contrast, the background thermal radiations do not pos¬ 
sess such features. 

Before ending this section, it is informative to comment 
on a key point to experimentally achieve the predicted 
microwave emission. In experiments, an Y 3 Fe 50 i 2 /Pt 
bilayer is conventionally grown on a millimeter-thick 
Gd 3 Ga 50 i 2 substrate. The thickness of the Y 3 Fe 50 i 2 
film is less than a micrometer in order for the effect of 
the spin transfer across the Y 3 Fe 50 i 2 /Pt interface to be 
visible. On the other hand, in our calculation we as¬ 
sume that, at low temperatures, the phonon mean free 
path grows to a millimeter scale (recall that in Ref. l3^ . 
from which we take the input data to draw Fig. [31 a bulk 
Y 3 Fe 5 0i2 was used). Therefore, an implicit assumption 
in our theory is that the interface between the Y 3 Fe 50 i 2 
film and the Gd 3 Ga 50 i 2 substrate is sufficiently trans¬ 
parent to phonons. This means that, although the qual¬ 
ity of Y 3 Fe 50 i 2 /Pt interface is usually the only matter of 
attention, much caution should be paid to the quality of 
Y 3 Fe 50 i 2 /Gd 3 Ga 50 i 2 interface to achieve the situation 
shown in Fig. [S] 


V. DISCUSSION AND CONCLUSION 

The physics behind the predicted microwave emission 
are best understood in terms of the spin-pumping pic¬ 
ture (see Fig. [I}. First, recall that the spin accumulation 
can be regarded as dynamically induced paramagnons in 
the language of many-body theory^^, and that the quan¬ 
tity Bqs. (I2|) and (jS]) denotes the den¬ 

sity of states of the paramagnons^!. Then, without the 
temperature bias, a magnon carrying spin angular mo¬ 
mentum —h pumps negative spin into P in the form of 
paramagnon excitations. This always results in an en¬ 
hancement of the Gilbert damping because the magnons 
lose the spin angular momenta [Fig. UKa)]. On the other 
hand, when P is heated, the pumped paramagnons tend 
to diffuse back into the cooler F, due to the drag effect of 
nonequilibrium phonons, which causes a forced spin back- 
flow against the pumped spin current [Fig. [IJb)]. This 
spin backflow gives a negative contribution to the Gilbert 
damping constant in accordance with Eq. (Ill|) under the 
condition that Qph = — Arph(Thot — Tcoid) <0. In this 
picture, the magnon instability in the limit of negligible 
intrinsic Gilbert damping constant a^°^ occurs when the 
sign of the spin-pumping current changes. 

Note that the present phenomenon differs from ther¬ 
mal spin-transfer torque in spin valve structures dis¬ 
cussed in Refs. [13 and jU as the spin valve structures 
are not required here. Furthermore, the phenomenon 
should be distin guis hed from a thermomagnonic analogue 
to Refs, m and in that the phonon-drag effect is the 
indispensable ingredient in the present theory. Moreover, 
the physics under discussion cannot be interpreted in 
terms of the spin-transfer torque exerted by a spin ac¬ 
cumulation caused by the spin Seebeck effect, because 
such a process would yield a result proportional to 
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(Jgjj required for accumulating spin and the other 
for exerting torque) instead of in the present theory. 
Note also that, since we do not consider spin valve struc¬ 
tures in this paper and there is only a uniform magneti¬ 
zation structure, we do not find any relevant processes in 
our calculation that are characterized by magnonic spin 
transfer torque. 

In conclusion, we have developed a theory for the ther¬ 
mal control of the Gilbert damping constant in a fer- 
romagnet/paramagnet bilayer, and predicted that a mi¬ 
crowave emission from the bilayer occurs at low tempera¬ 
ture, where the participation of nonequilibrium phonons 
is remarkably enhanced. Since the structure of the bi¬ 
layer is much simpler than that of spin valves, the present 
phenomenon extends the range of applications of spin 
caloritronics. 
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Appendix A: Derivation of Eq. (1111) 


interaction across the interface, we obtain 


Qph — 


hL 
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where L = (fl(i’ 2 )^( 2 . 3 )) 2 ^.a. 2 )^G, 3 ) 
the number of lattice sites at the di/dj interface, and 
we introduced the shorthand notation f = We 

evaluate the integral over ly by picking up the dominant 
contribution that is proportional to the phonon lifetime 
Tph, and obtain 


Qph 


fiLxph ^ r'K2 Iin^Ki(r'K2)Iin-DK3(;^K2) 
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where the result was linearized with respect to AT. After 
a little algebra, Qph is calculated to be 


Qph 

i^ph 


-ifphAT, 

7r%_ Pph(’^p){Pph(’^p)y{kBT)^ 
{hvpfT 
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where vp are the Debye frequency of P, Pp^(p) (Pph(p)) 
is the density of states of phonons in P (F), and the 
constant cg is defined by 


In this section, we derive Eq. dill) in the main text. In 
the following, it is covenient to introduce a notation di 
{i = 1,2, 3) to represent a temperature domain in Fig. [2] 
having its local temperature T^. 

We begin by deriving the expression of the phonon heat 
current Qph flowing into P. This quantity can be calcu¬ 
lated in a manner similar to Ref. using the definition 
Qph = (^-ffph^); where (• • •) means the statistical aver¬ 
age, and TJph^ = ^K 3 Sk 3 ^K 3 Hamiltonian 

for phonons in domain dg with 6 x 3 being the phonon an¬ 
nihilation operator defined in the domain. To describe 
heat transfer across the domains, we need to introduce 
the interaction Hamiltonian 

^?h'^ = ^ E +H.c.(Al) 

K,.K, 

where Hk* = ^Ki + is the Fourier 

transform of ErpGinterface - i-o) with 

being the characteristic energy of the phonon transfer 
across the di/dj interface. 

(3) 

Using the Heisenberg equation of motion for and 
then performing perturbation approach in terms of the 
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Appendix B: Derivation of Eq. dl3l) 


In this section, we derive Eq. (US in the main text. 
The spin current injected into P via the spin Seebeck 
effect is given by^ where 

s{ri,t) = Sk(t)e‘'^''’‘ is the spin accumulation 

in P. Using the Heisenberg Eq. of motion for and as¬ 
suming the steady state condition, this quantity can be 
calculated from = X^k.q L Re[Jsd(k - q)Gi^q(u;)], 
where Gj^q(a;) is the Fourier transform of the interfacial 
Green’s function = —i(aq(t')sk(t)), and we in¬ 

troduced the shorthand notation f = 

Now we focus on the process shown in Fig. [IJd). Fol¬ 
lowing the procedure given in Ref. l4ll the corresponding 
spin current injection is represented as 


rSSE 



-Fk,q(j^)ImJT>^3(p), 
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where 5D^^{v) was defined in Eq. ([7]), and 
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with k_ being introduced below Eq. ([5]). We 

first use the relation coth(A — B) — coth(A) = 
sinh(i3)/[sinh(A) sinh(A — B)] and then perform the in¬ 
tegral over uj using ImGq(a;) = —'!rS{uj — Wq). Then, 
is calculated to be 


2TTkBTi^Ts{ C‘'2k^T'’) sinh(^^) 


-Fk.q(i^) = - 




sinh sinh I 

ImXk(w) 




hiO, 
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where the derivative with respect to w originates from 
the fact that only the even-in-i/ component gives a non¬ 
vanishing contribution to because of the symmetry 
of the V integral. Substituting Eq. (IB3I1 into Eq. (jBip and 
performing the integral over v by picking up the phonon 
poles, we obtain Eq. m, where we have used an approxi¬ 
mation sinh(~ 1 because the dom¬ 
inant contribution comes from a region ^ 


Appendix C: Comparison of the theory with 
experiments 

Equation m means that the Gilbert damping con¬ 
stant is reduced when the paramagnet side is hotter, the 


sign of which is consistent with the experiment reported 
in Ref. Q For its order of magnitude estimation, we use 
Eq. (fTHll . Since we have ~ 2.0 x 10"^ at room tem¬ 

perature Tr under a temperature bias of AT = 20 K^, 
we have an estimate Sa' AT/[a^°'> + 6a] ~ 30%, where we 
used the fact that <C 6a. This estimate is compara¬ 
ble to the experimentally observed reduction of 17%. 


On the other hand, one could attribute the observed 
reduction to a trivial effect of the temperature bias be¬ 
cause there might be a temperature bias effect that is 
free from the phonon heat current Qph and comes solely 
from temperature dependence of material parameters in 
Eq. ([2|). To see this, we first use the fact that the dom¬ 
inant temperature dependence comes from that of the 
paramagnetic spin susceptibility and the size of the 
spin (oc saturation magnetization where the for¬ 

mer quantity is defined as a function of the tempera¬ 
ture of P (i.e., T = Tp) and the latter as a function 
of the temperature of F (i.e., T = Tp). We next use 
the relation Tp = Tb + \A.T and Tp = Tb — ^AT 
by setting Tr = ^{Tp + Tp), and calculate (5a (Tr -|- 
AT) = const X Ms(Tr -|- ^AT) x xp{Tb — ^^T). Then, 
we obtain [i5a(TR -|- AT) — (5a(TR)]/[a^°^ -I- 5a(TR)] m 
i^[lnxp(T)-lnM,(T)]T =TrAT, where we have again 
used the fact that a^^^ <C (5a(TR). This estimate gives 
at most a 2% reduction in a under a temperature bias of 
AT = 20 K, which is much smaller than the experimen¬ 
tally observed value. This means that Eq. (HU is indeed 
the origin of the observed behavior. 
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